Calculation of Excavation Force for ISRU on Lunar Surface 
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ABSTRACT 

Accurately predicting the excavation force that will be encountered by digging tools on the lunar surface is a 
crucial element of in-situ resource utilization (ISRU). Based on principles of soil mechanics, this paper develops an 
analytical model that is relatively simple to apply and uses soil parameters that can be determined by traditional soil 
strength tests. The influence of important parameters on the excavation force is investigated. The results are 
compared with that predicted by other available theories. Results of preliminary soil tests on lunar stimulant are also 
reported. 

Nomenclature 

a h = acceleration in horizontal direction (m/s 2 ) 

a v = acceleration in vertical direction (m/s 2 ) 

c = soil cohesion (N/m 2 ) 

c a = soil-tool adhesion (N/m 2 ) 

d = depth of excavation (m) 

F side = side friction (N) 

Fbiade = friction on the blade (N) 

g = gravitational acceleration (m/s 2 ) 

G s = specific gravity 

H r = corrected depth of fall (cm) 

K = Gill’s cut resistance index (N/m) 

K pe = dynamic passive earth pressure coefficient 
K 0 : = earth pressure coefficient at rest 

L w = length of the failure wedge at surface (m) 

P p = passive earth pressure (N) 

q = surface surcharge (N/m 2 ) 

T = total excavation force (N) 
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T = temperature (deg) 

t = time (minute) 

v = tool velocity (m/s) 

w = width of excavation blade (m) 

W = tool length (m) 

W b = weight of excavation blade (N) 

W s = weight of soil wedge (N) 

a = inclination angle of the blade (deg) 

Op = inclination angle of the failure wedge (deg) 

(3 = inclination angle of the side friction (deg) 

pi = external friction angle (deg) 

y = unit weight of soil (N/m 3 ) 

6 = friction angle between soil and blade (deg) 

c|) = internal friction angle 

cp = inclination angle of acceleration (deg) 

p = mass density of soil (kg/m’) 

p w = mass density of water (g/cm 3 ) 

r| = viscosity of water (g/cm-s) 

9 = shear plain failure angle (deg) 

a v = vertical normal stress in soil (N/m 2 ) 
cth = horizontal normal stress in soil (N/m 2 ) 

I. Introduction 

O ne of the crucial operations for the next manned mission to the moon is the in-situ resource utilization (ISRU), 
which requires to obtain critical energy sources and life supporting materials such as oxygen and hydrogen 
from the surface regolith on the moon. In order to design excavation tools for digging and pushing lunar regolith for 
in-situ resource utilization (ISRU) applications, the excavation force that may encountered by the tools need to be 
accurately calculated. Over the years, four two-dimensional models have been developed to estimate the excavation 
force. These are the ones developed by Gill and Vanden Berg 1 , Swick and Perampral 2 , McKyes and Edwards 3 , and 
the Martin-Marietta Corporation 4 in association with the Colorado School of Mines for the design of the Viking 
Mars exploration probes. 

Even though these models have been developed and applied quite frequently, most of them are empirical and use 
parameters difficulty to determine in laboratory. Therefore, there is the urgent need to develop a model to calculate 
the excavation force based on principles of soil mechanics, which should be relatively easy to use and parameters be 
easy to obtain using conventional soil tests. In addition, this model needs to be validated by other commonly used 
engineering methods such as finite element method and laboratory measurements. In this paper, an analytical model 
that is relatively simple to apply and uses soil parameters that can be determined by traditional soil strength tests is 
developed. The influence of important parameters on the excavation force is investigated. The results are compared 
with that predicted by other four available models. 

II. Results of Preliminary Tests on Simulants 

Since only a few kilograms of lunar soil were brought back from the moon, simulants that have properties similar to 
that of the lunar regolith need to be created so that engineering tests can be conducted to study the conditions that 
would be expected on the moon. Currently the most commonly used simulant is the ISCla, which has three sub- 
groups: ISCla fines, ISCla very fines, and ISCla. Particle size analysis was conducted on the fines in ISCla fines 
and very fines at the geotechnical laboratory of Case Western Reserve University. Hydrometer test was used to 
determine the grain size distribution of fine-grained particles with sizes smaller than 75 pm and down to 1pm. The 
principle of hydrometer test is based on Stokes’ law, which gives the particle diameter D (mm) as 

D = V{[30t 1 H r ]/[ 981(G s - l)p w t]} (1) 

in which t is the time (min) after the beginning of sedimentation, G s the specific gravity of the soil particles, p w the 
unit mass of water (g/cm 3 ) at temperature T, r| the viscosity of water (g/cm-s) at temperature T, and H R the corrected 
depth of fall (cm). Since most soil particles have irregular shape and EQ. 1 was derived for spheres, the particle size 
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determined from this method is the equivalent spherical diameter of the particles. The test has been standardized by 
the American Society for Testing and Materials (ASTM) as ASTM D422. A test takes about 3-5 days to finish. The 
result on JSC la very fines (VF) is shown in Figure 1. 

As shown in the figure, all the fine particles have particle sizes smaller than 25 pm. There is a wide range of 
particle size with D 10 = 1.2 pm, D 30 = 3 pm, and D 60 = 6pm, where D 10 , D 30 , and D 60 means that there are 10%, 30%, 
and 60% (by weight) of the particles have a size smaller than the number, respectively. The coefficient of uniformity 
= D 60 /Dio = 5 and coefficient of curvature = D 3 o 2 /(D 10 -D 6 o) = 1-25. There are about 6% of the particles with a size 
smaller than 1pm. Based on particle size analysis, the fines have about 50% silt and 50% clay. The soil is classified 
as a well-graded clay. 






Hydrometer Analysis 
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Figure 1. Particle size distribution of JSCla VF 


For the JSCla fines (F), one sample each was taken at the top and bottom of the storage bucket, respectively. 
The sample at the top was spooned out the bucket while the sample near the bottom was taken by a mechanical thief. 
The objective is to see whether there is segregation of particles during the transportation of the fines. The results of 
the hydrometer tests on JSCla fines are shown in Figure 2. 

As shown in the figure, the particle size distribution of the two samples are almost identical, indicating there is 
no segregation of particles during the transportation process. There is a wide range of particle size with Di 0 = 6.2 pm, 
D 30 = 12 pm, and D 60 = 2 1pm, where D 10 , D 30 , and D 60 means that there are 10%, 30%, and 60% (by weight) of the 
particles have a size smaller than the number, respectively. The coefficient of uniformity = D 60 /Di 0 = 3.4 and 
coefficient of curvature = D 30 2 /(D 10 -D 60 ) =1.1. There are about 2% of the particles with a size smaller than 1pm. 
Based on particle size analysis, the fines are classified as poorly graded silt. 

The specific gravity of the fines is determined using the method described by ASTM D854. It involves the use of 
Archimedes’ principle to determine the volume of irregularly shaped soil particles. The specific gravity is the ratio 
of mass density of the solid particles divided by the mass density of water at 4°C. The results of specific gravity test 
show that the specific gravity of the JSCla fines is 2.85. This is typical for silt and clay particles. 
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Figure 2. Particle size distribution of JSCla F 

TIT. A Model to Calculate Excavation Force Based on Principles of Soil Mechanics 

This model to calculate excavation force on an excavation blade is based on the principles of soil mechanics. The 
problem is illustrated in Figure 3. 

surcharge q 



Figure 3. Excavation blade and soil body at failure 

The free-body diagram of the blade is shown in Figure 4. 
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Figure 4. Forces act on the blade 

Assuming the blade and the body of soil in failure state have a horizontal acceleration of a h and a vertical 


acceleration of a v , the equations of motion for the blade can be written as: 

T x + F b i ad e sina - P p cos (a-8) - F side cosp = (W b /g)a h (2) 

T y - F blade cosa - W b - P p sin (a - 8) - F side sinp = (W b /g)a v (3) 

in which T x and T y are the horizontal and vertical components of T, respectively, and p is the inclination angle of the 
side friction force. The equations can be re-written to calculate T x and T y : 

T x = -F b i ade sina + P p cos (a-8) + F slde cosP + (W b /g)a h (4) 

T y = F blade cosa + W b + P p sin (a - 8) + F side sinp + (W b /g)a v (5) 

Then, the total excavation force T is 

T = V(T X 2 + T y 2 ) (6) 


Three forces need to be calculated in these equations: friction force on the blade F blade , passive earth pressure P p , 
and side friction force F side . In general, since the blade is moving away from the soil behind the blade, the normal 
stress on the surface of the blade is expected to be small, resulting in a small frictional resistance. For soil with small 
cohesion, such as the surface soil expected on the moon, the cohesion component of the friction is also small. 
Therefore, F b i ad e can be ignored without much influence on the calculation. 

For calculation of P p , the widely used Mononobe-Okabe 5 theory to calculate dynamic earth pressure on a 
retaining wall during earthquake is modified for this application. It is a pseudo-static type of calculation in which the 
dynamic effect is taken into account by inertial forces. For the soil wedge in front of the blade shown in Figure 5, 


assuming a failure state, the passive earth pressure P p on the blade is given by 

P p = 0.5 K PE (1 + a v /g) yd 2 W + 2cdwVK PE + K PE qdW (7) 

in which K PE is the dynamic passive earth pressure coefficient which is given by 

cos 2 (c|) + a + <p) 

K pe = ----- - - - (8) 


cos<pcos 2 acos(S- a - q>){ 1 - V[sin(8 + cj))sin((j) +(p)]/[cos(S - a - <p)cosa]} 2 
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where the angle (p is the inclination angle of the total acceleration 

cp = tan _1 [a h /(g+a v )] (9) 

The critical failure surface is inclined from horizontal by an angle 

oc p = -cp - (> + tan 1 {[tan(<t> - cp) + C 3E ]/C 4E } (10) 

where 

C 3E = V{tan((|)+(p)[tan((|)+(p) +cot(c|)+a+(p)][l+tan(5-(p-a)cot((j)+a+(p)]} (1 1) 

C 4E = 1 + {tan(S-cp-a)[tan(<|)+(p) + cot(c|)+a+(p)]} (12) 


For the side friction calculation, assuming the failure area is the same as the failure wedge in front of the blade and 
the direction of movement the same as the inclination angle of the wedge or (5 = a p , then as shown in Figure 6, the 


vertical normal stress at depth y is 

G v = Yy + q (13) 

The horizontal normal stress at the same level is 

G h = K 0 g v = K 0 (yy + q) (1 4) 

The length of the failure wedge at the surface is 

L w = d (tana + cotap) (15) 

At depth y, the length of the failure wedge is 

L(y) = L w [(d-y)/d] (16) 
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surcharge q 



Figure 6. Side frictional force on the soil wedge 

For a small slice of the wedge at depth y with a thickness of dy, the side friction force is 

dF slde = 2(c + o h tan(|))L(y)dy = 2 (c + K 0 yy tancj) + K 0 q tan<|)) L w [(d-y)/d]dy (17) 

For the entire side surface 
d 

F s ide = 1 2 (c + K 0 yy tan<|) +K 0 q tan<» L w [(d-y)/d]dy 
0 

= L w (cd + K 0 qd tancj) +K 0 y tanc[)d 2 /3) (18) 

Substituting P p and F side back into equations (4), (5), and (6), the excavation force T can be determined. 

IV. Parametric Study and Interpretation 

The model described above (herein called Zeng model) and the other four existing models can be easily 
programmed using Excel and the influence of important parameters be studied. Parametric studies were done to 
demonstrate how these five models respond to variances in their inputs. This gives insight into what are the most 
important soil parameters for design. The parametric studies hold all inputs constant except for one. The variable 
input is graphed along the x-axis with the draw bar force along the y-axis. The actual base values for all these inputs 
are the same as those used in the 2006 Wilkinson-DeGennaro 4 paper on regolith excavation and traction. Most of 
these values were taken from the Lunar Sourcebook and represent feasible values for the lunar soil. 


Tool Width (w): 1 m 

Tool Length (W): 0.7 m 

Tool Depth (d): 0.5 m 

Tool Velocity (v): 0.1 m/s 

Moon Gravity (g): 1 .63 m/s 2 

Horizontal Acceleration (a h ) 0 m/s 2 

Vertical Acceleration (a v ): 0 m/s 2 

Soil Surcharge (q): 1 kg/m 2 
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Soil Density (p): 1680 kg/m 3 

Gill’s Cut Resistance Index (K): 1000 N/m 
Soil Cohesion (c): 170 N/m 2 

Soil-Tool Adhesion (c a ) 1930 N/m 2 

Rank Angle (a) 45 deg 

Internal Friction Angle (<])) 35 deg 

External Friction Angle (Pi) 10 deg 
Shear Plain Failure Angle (9) 30 deg 

Soil-Blade Friction Angle (5) 20 deg 

At Rest Earth Coefficient (K 0 ) 0.573 

There are several characteristics that are unique to the Zeng model. Most notably is its ability to handle 
acceleration of the tool blade while the other models assume the soil mass is moving at a constant speed. The Zeng 
model is also the only model to calculate passive earth pressure. This passive earth pressure must be inclined at a 
certain angle, which is prescribed. It is often between one half and two-thirds the value of the internal friction angle. 
The input of this soil-blade friction angle for determining the line of action of passive earth pressure is required. 
Parametric studies for the Zeng model, broken into its passive pressure and side friction components are shown in 
Figure 7. It shows how horizontal acceleration, vertical acceleration, the soil-blade friction angle, and external 
friction angle affects the total excavation force, the passive earth pressure components, and the side friction 
component, respectively. 

Parametric studies of the five soil models are conducted concurrently and the results are shown in Figure 8. 
Note that the Swick and Perumpral and McKyes models are quantitatively equal and hence will be represented by 
the same line throughout this paper. Regarding the graph in the upper left, the Zeng model does not consider the 
velocity of the toil/soil mass and is thus a horizontal line, while for the other models the value of the draw bar force 
varies, the way the draw bar force changes is similar. It increases deliberately with increasing velocity, at least on 
the range shown in the graph. Even at very low speeds, the draw bar force does not differ significantly from its 
value at moderate speeds. The graph in the upper right, a parametric study of force versus gravity, shows how the 
models would work under different planetary gravitational fields. Similar to velocity graph, there is a dichotomy of 
the models, with Zeng and Swick/McKyes models in one group and the Gill and Viking models in another. Another 
similarity to the velocity; while the models are separated into two sets with similar nominal values, the change in 
force is a general increase in both cases. 

The next two graphs are parametric studies of variance with tool depth and soil density. Both these 
parameterizations also show the dichotomy between the Zeng/Swick/McKyes models and the Gill/Viking models. 
Also of note in these two studies, as well as the gravity graph, are the following characteristics. The Zeng model 
starts out at a lesser value than the Swick/McKyes but increases at a slightly faster rate. Similarly the Viking model 
starts out at a smaller value than the Gill model but increases moderately quicker. Also, not only do the Gill/Viking 
models start out at higher nominal values than the Zeng/Swick/McKyes models, but the differences between them 
increase with the increasing parameter. Variances in soil density and gravity both show linear changes in draw bar 
force, while variances in soil depth lead to quadratic or quasi-quadratic increases in draw bar force. 

The last two charts in Figure 8 show parameterizations of the internal and external friction angles. The Viking 
model does not consider either friction angle, so it remains a straight line for in both graphs. For changes in both 
angles, the Gill model gives factor of 2 like differences for the draw bar force as compared to the Swick/McKyes 
and Zeng models. The Gill model also shows the previously identified pattern of increasing at a greater rate. The 
Swick/McKyes models behave in a similar fashion to what they did in the previous models but the Zeng model 
responds differently to angle changes. With changes in internal friction angle the Zeng model starts out at a similar 
place to the Swick/McKyes models at low friction angles, but increases slower. In the external friction angle 
parametric study the Zeng model starts out at a similar place and decreases slightly in value as the angle increases. 

These last two parametric studies, shown in Figure 9, are the most unique, particularly in the behavior of the 
Zeng model. The other four models exhibit behavior similar to what we have seen before. The Viking and Gill 
models give higher values for horizontal force than the Swick/McKyes models. The one notable break from that 
standard is the linear behavior of the Viking model’s variability with rank angle. The Zeng model’s dependence on 
cohesion is much stronger than any of the other models. It shows a rapid increase in the required drawbar force as 
cohesion increases, having the lowest values for cohesionless soil but the highest values for very cohesive soil. This 
strong dependence on cohesion shown by the Zeng model might help describe some of the unusual results seen in 
lunar excavation. Assignment of a value for the apparent lunar cohesion still remains a delicate process. The Zeng 
model’s behavior for rank angle is stranger still, showing a general decrease in force as rank angle increases at small 
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angles. As the angle increases, the slope becomes flatter, reaching at minimum at a rank angle of about 57 degrees. 
Beyond this minimum, horizontal force begins to increase again. 
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Figure 7. Influence of important parameters on excavation force 
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Horizontal Force vs. Velocity 


Horizontal Force vs. Gravity 
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0.4 0.6 

Depth (m) 



Soil Density (kg/m 3 ) 



Figure 8. Comparison of results of excavation forces by four different models 
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Horizontal Force vs. Cohesion 
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Figure 9. Influence of cohesion and rake angle on excavation force 


V. Conclusions 

Based on the results of this study, the following conclusions can be drawn: 

1) The particle size analysis has been successfully carried out on JSC la simulant. The fines and very fines of 
the simulant has been characterized 

2) A new model to calculate the excavation force has been developed. This model is based on the principles 
of soil mechanics and quite easy to apply. The parameters used in the model can be determined by 
conventional soil tests. Important parameters that influence the magnitude of the force are identified. 

3) Parametric study shows that excavation force calculated by this model is more sensitive to cohesion of the 
regolith, which may help to explain the high value of excavation forced measured in past mission to the 
moon. 
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